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We consider v^, — > v x oscillations in the context of the Mass Varying Neutrino (MaVaN) model, where the 
neutrino mass can vary depending on the electron density along the flight path of the neutrino. Our analy- 
sis assumes a mechanism with dependence only upon the electron density, hence ordinary matter density, of 
the medium through which the neutrino travels. Fully-contained, partially-contained and upward-going muon 
atmospheric neutrino data from the Super-Kamiokande detector, taken from the entire SK-I period of 1489 
live days, are compared to MaVaN model predictions. We find that, for the case of 2-flavor oscillations, and 
for the specific models tested, oscillation independent of electron density is favored over density dependence. 
Assuming maximal mixing, the best-fit case and the density-independent case do not differ significantly. 
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PACS numbers: PACS numbers: 14.60.Pq, 14.60.St, 96.50.Sf 
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I. INTRODUCTION 

Neutrino oscillations result when at least one neutrino has 
mass different from the others, so that neutrinos mass eigen- 
states are distinct from the flavor eigenstates. For atmospheric 
neutrinos, — > V x oscillations f 1] are strongly favored over 
Vjj — » v. v oscillations |2] and other exotic mechanisms for neu- 
trino disappearance y|]. Neutrino oscillation without the in- 
troduction of sterile flavors has also been sufficient to resolve 
the solar neutrino problem 0, \^ . 

Super-Kamiokande previously reported best fit parameters 
for atmospheric two-flavor — ► V x oscillations, providing 
an explanation of the atmospheric neutrino anomaly 16J]. In 
Super-Kamiokande, where previous analyses always consid- 
ered only geometric path length independent of medium, neu- 
trino disappearance effects are strikingly evident for upward 
going neutrinos, which have passed through many km of 
rock, in comparison with downward-going neutrinos, which 
travel mainly through air. We note that other experiments 
that observed neutrino oscillations, such as KAMLAND |7], 
K2K H, and MINOS [9], detected neutrinos whose path was 
almost entirely through rock. 

Here we consider a possible consequence of the Mass Vary- 
ing Neutrino (MaVaN) model ifioll . In this model, the neutrino 
mass can vary depending on the matter density along the path 
of the neutrino. MaVaNs could provide a source for the dark 
energy fi"Hl . Here we will assume that the mass variation of 
the neutrinos depends only upon the electron density of the 
environment, a possible side effect from radiative couplings 
of active neutrinos and electrons II 1 Oil . (In ordinary matter 
comprising the Earth's atmosphere and interior, overall matter 
density and baryon density, considered in other MaVaN sce- 
narios, are essentially proportional to electron density.) If this 
hypothesis is accurate, it may be possible to probe the mass 
dependence with Super-Kamiokande data. 

We used atmospheric neutrino data from the Super- 
Kamiokande-I running period (1996-2001) to test whether 
V,u — > v x oscillations have an apparent dependence on the 
electron density of the material the neutrino passes through. 
For this analysis, we assumed that the mass squared differ- 
ence is proportional to some power of the electron density, 
Amj^y ~ p" . We also assumed that the mixing angle is con- 
stant for all media. 



H. DATA ANALYSIS 

Super-Kamiokande is a water Cherenkov experiment lo- 
cated within Mt. Ikeno-yama in central Japan, under 2700 
meters water equivalent rock overburden. It has a cylindri- 
cal design, holds 50 kilotons of water, and is divided into 
two optically separated sections by a structural framework 
that supports the photomultiplier tubes (PMTs). During the 
SK-I running period from which the data analyzed here came, 



the detector had an inner detector (ID) equipped with 1 1 146 
50cm PMTs aimed inward, and an outer detector (OD) vol- 
ume instrumented with 1885 20cm PMTs, aimed outward and 
equipped with wavelength-shifting plastic plates. The OD 
functions primarily as a veto counter, tagging charged parti- 
cles that enter or exit the ID. Within the ID we define a cen- 
tral 22.5 kiloton fiducial volume, within which detector re- 
sponse is expected to be uniform. Fully-contained (FC) neu- 
trino events are those where interaction products are observed 
in the ID, with no significant correlated activity in the OD, 
while partially-contained (PC) events are those where some 
interaction products exit the ID. Upward-going muon events 
are those where a penetrating particle travelling in the upward 
direction enters and either stops or passes through the detector, 
and are attributed to muons produced by neutrino interactions 
in the surrounding rock. In general terms, FC, PC, and up- 
ward muon events represent successively higher energy sam- 
ples of neutrino interactions, ranging from 200 MeV for the 
lowest energy FC events to above 1 TeV for the highest energy 
upward-going muons. Further details regarding the Super- 
Kamiokande detector design, operation, calibrations, and data 
reduction can be found in Il6l [l2ll . 

Super-Kamiokande-I (SK-I) data taking for physics analy- 
sis began on May 17, 1996 and lasted until a planned shut- 
down for refurbishing on July 16, 2001, including a total of 
1489.2 days of effective livetime (92 kiloton-years) for FC and 
PC events, and 1645.9 days for upward-going muon events. 
The full Monte Carlo (MC) event sample generated for com- 
parison is equivalent to a 100 live-year period. The SK-I 
database used in the analyses presented here has the following 
statistics: for FC events, 12180 data and 13676.7 MC events in 
a livetime-scaled sample, for PC events, 911 data and 1 129.6 
MC events. The upward-going muon statistics are: for stop- 
ping muon tracks (those that enter but do not exit the detector), 
417.7 data and 713.5 MC events, and for through-going muon 
tracks, 1841.6 data and 1669.5 MC events, after background 
subtraction for near-horizontal downward-going muons as de- 
scribed in ref . |0] . 

Two different analyses are considered here. The first anal- 
ysis does not take into account path length in air for neu- 
trino oscillations. Instead, it assumes MaVan-type oscillations 
only occur in high-density matter, taking into account surface, 
crust, mantle, and Earth core densities. In other words, the 
neutrino flight path length is in effect taken to be only that 
portion of the geometric path that lies in high density mat- 
ter. Because the path length systematic error is not relevant 
in this test, it is removed from the list of 39 systematic errors 
that was used in previous Super-Kamiokande analyses. The 
second analysis also takes into account path length in air and 
its density in considering neutrino oscillations effects in the 
context of the MaVaN model. 

A detailed description of the SK-I atmospheric neutrino 
chi-squared zenith angle analysis using the "pull" method lfl3ll 
can be found in [6], and the same methods are applied here. A 
profile of the mountain surrounding Super-Kamiokande was 
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taken from topological maps to determine the downward neu- 
trino path length in rock. The maps used are the United States 
Geological Survey agency's Digital Elevation Maps lfl4ll . with 
data points at 7.5 minute (about 30 meters) spacing. 

In the conventional 2-flavor — * V x oscillation framework, 
the oscillation probability can be written as 



sin 2 29 sin 2 



1.27Am 2 L GeV 



km-eV" 



where L is the distance travelled by the neutrino between pro- 
duction and detection, and E v is the neutrino energy, while the 
oscillation parameters are the mixing angle 9, and the mass 
difference squared, Am 2 = \m^_ — m\ | in the usual nomencla- 
ture. To test for density dependence, we replaced Am 2 in 
the equation above with an effective mass difference that is 
proportional to the electron density of the medium, Am 2 — ► 

Am 2 x (jf*j > where p e is the electron density of the matter 

in neutrino trajectory, and n parameterizes the density depen- 
dence. For neutrinos passing through layers of matter with 
different density, the path length in each layer is taken into 
account, po is set at 6.02 x 10 23 e/cm 3 . Several different val- 
ues of n are tested here, as described below. To approximate 
the mass density of the Earth for upward traveling neutrinos, 
the Preliminary Earth Reference Model (PREM) 1 15] is used. 
The electron density is taken to be the mass density multiplied 
by the charge-to-mass ratio lfl6ll . 

Atmospheric neutrinos are produced at altitudes around 
15 km above sea level. A neutrino travelling nearly verti- 
cally downward thus passes through about 15 km of air (low- 
density matter) followed by 1-2 km of rock (high density 
matter) before interacting within the detector, while a nearly- 
horizontal downward-going neutrino travels through over 200 
km of air and about 10 km of rock. Thus, a downward neu- 
trino has 8-20 times as much path length in air as in rock. 
If oscillation only occurs in high density matter, the effective 
path length will be shorter than the geometric path length trav- 
elled by the neutrino from point of production to detection. In 
that case, the standard oscillation model, which uses the full 
geometric path length, assumes more oscillation cycles occur 
for a given Am 2 , and oscillation effects are more significant 
for neutrinos at low energies. Analyses both including, and 
neglecting, the air path length are considered here. 



A. Analysis Neglecting Air path length 

Density dependent oscillation models were first tested ne- 
glecting the air path length and using specified values of n 
suggested by theorists 11711 . The results are shown in Table U 
for two cases: fits constrained to lie in the physical region, 
where maximal mixing is assumed for all densities, and with 
this constraint relaxed, where mixing angle is allowed to vary. 
All models tested are ruled out, with confidence level better 
than 99.9%, when compared to fits assuming oscillations in- 
dependent of density of the medium (including air). 

Next, we treat the density dependence as a free parameter, 
and find the best fit electron density dependence exponent n 



for the case of matter-dependent oscillations assuming maxi- 
mal mixing (sin 2 (29) = 1) in all densities (except air). This 
analysis uses both n and Am 2 as the varied parameters. Am 2 
values from 10~ 4 to 10° eV 2 were used over 81 bins, and den- 
sity dependence powers from n — — 3 to +3 were considered, 
over 61 bins. The results produce a y^ njn = 187.3/178 d.o.f. 
for (n,Am 2 ) = (-0.30,3.16 x 10~ 3 eV 2 ). Thus, the model in 
which air path length is neglected, i.e., oscillations occur only 
in rock, is disfavored at about the 3.5 O level relative to stan- 
dard oscillations. 



B. Analysis Including Air path length 

We repeated the analysis, this time taking into account the 
path length in air. Now the n = case, where there is no de- 
pendence on the electron density, corresponds to the conven- 
tional oscillations model, since the full geometric path length 
is used. 

The air mass density depends greatly on the altitude of the 
air above sea level. For simplicity, a constant value of 10~ 3 
g/cm 3 (corresponding to actual air density at mountain alti- 
tude) was used for the mass density. 

The same set of fixed n values are tested again, to check 
various density dependences, this time including the air path 
length. The results are listed in Table [II] All models 
tested are excluded by > 3.8a when compared to the density- 
independent case. 

Finally, we again treat density dependence n as a free pa- 
rameter, this time taking air path length into account, once 
more assuming maximal mixing. The results are shown in 
Figure Q] A minimum % 2 = 174.3/178 d.o.f. is found at 
(n,Am 2 ) = (-0.04, 1.95 x 10~ 3 eV 2 ). This result and the stan- 
dard 2-flavor oscillation result (i.e. with n = 0) do not differ 
significantly. 



III. SUMMARY AND CONCLUSIONS 

We used data from Super-Kamiokande-I to test mass- 
varying neutrino models, by considering oscillation effects 
where the effective mass of a neutrino is dependent on the 
matter environment of the neutrino flight path. Our analy- 
sis assumed varying neutrino mass was due entirely to the 
electron density of the material in a neutrino's path, a pos- 
sible consequence of the MaVaN theory. Assuming a constant 
mixing angle in all environments, we substituted Am 2 with 

Am 2 (jf^j (where p„ = 6.02 x 10 23 e/cm 3 ) for specific val- 
ues of n suggested by theorists, as well as with n as a fitted 
parameter. Two different models were considered. In the first 
analysis, we neglected the portion of the neutrino flight path 
in low density matter (air), and assumed mass-varying effects 
occur only in high density matter (rock). For the second anal- 
ysis, we allowed oscillation effects in all portions of the neu- 
trino path, but (in contrast to the conventional oscillation anal- 
ysis, where only the geometric path length is considered) took 
into account mass densities. The primary difference in the two 
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TABLE I: Comparisons of % 2 Values for different models without air path. There are a total of 178 degrees of freedom in the % 2 (see ref. ((J). 
Note that n = is not equivalent to oscillations independent of medium, because here we do not allow for any oscillations in air. 
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TABLE II: Comparisons of % 2 Values for Different Models including air path. Note that here, n = is equivalent to oscillations independent 
of medium. 
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hypothesis lies in the effective path lengths for near-horizontal 
downward-going neutrinos. 

In the "pull" fitting procedure used, some of the systematic 
errors are allowed to vary in the fit. When fitted systematics 
are forced significantly outside a reasonable range of values, 
it is an indication that the model used is disfavored by the 
data. For the high-density-only oscillation model, some of 
the systematics, particularly the V^/v e ratio for energies be- 
low 5 GeV, adjusted to values beyond a reasonable deviation 
from the expected value. This does not happen when fitting 
to the conventional oscillation model, and results in a larger 
chi-squared value for the mass-dependent models. Based on 
relative chi-squared values in our results, the hypothesis that 
neutrinos oscillate only in high density matter is disfavored 
relative to conventional oscillations by at least 3.5 a, without 
requiring additional effects. 

For the second analysis, where oscillations in low density 
matter were included, a mass difference with electron den- 
sity dependence is disfavored at more than the 3.8-a level 
when compared with neutrino oscillations with a fixed mass 
difference, for all fixed values of n tested. In addition, a 
freely varying density dependence analysis was performed, 
assuming maximal mixing. It produced a minimum chi- 



squared value of xl in = 174.3/178 d.o.f. at (« = -0.04, Am 2 = 
1.95 x 10~ 3 eV 2 ), consistent (within 1.4 o) with the n = (no 
density dependence) case. 

The results of our analyses show no evidence that the en- 
vironmental electron density influences the effective Am 2 de- 
termined using Super-Kamiokande atmospheric neutrino data. 
We have not explicitly considered models where the mixing 
angle is not constant in all densities, or which assume 3 flavor 
oscillations, nor do we exclude all variations of MaVaN mod- 
els. We find that conventional density independent — > V x 
oscillations are sufficient to explain the atmospheric neutrino 
data. 
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FIG. 1: Am 2 — » Am 2 x (j^^J (including air path length). Upper 

plot shows relative-/ 2 confidence level contours on the Am 2 versus n 
plane, obtained when taking into account both high and low density 
matter path lengths. The lower plots display the y} — % 2 nin contours, 
with confidence levels shown, at the best-fit parameter values. 



